in 
o 
o 

(N 
S3 . 

Possible crossover of a non universal quantity 
cn : at the upper critical dimension 

S. Galam 1,2 and A. Mauger 1 ' 3 
^aboratoire des Milieux Desordonnes et Heterogenes 
Case 86, 4 place Jussieu, F-75252 Paris Cedex 05 
Centre de Recherche en Epistemologie Appliquee, 
CREA - Ecole Polytechnique, CNRS UMR 7656, 
1, rue Descartes, 75005 Paris 
t3 ! 3 Federation de Recherche Matiere et Systemes Complexes, 

CNRS UMR, 140, rue de Lourmel, Paris 75015 



CO 

c3 



> 



C3 



[galam@ccr.jussieu.fr; mauger@ccr.jussieu.fr 



Abstract 



We report on a possible crossover of a non universal quantity at 
ly-j ■ the upper critical dimensionality in the field of percolation. Plot- 

■ ting recent estimates for site percolation thresholds of hypercubes in 

dimension 6 < d < 13 against corresponding predictions from the 
GM formula p c = po[(d — l)(q — \)}~ a d b for percolation thresholds, a 
significant departure of p c is observed for d > 6. This result is rem- 
iniscent of the crossover undergone by universal quantities in critical 
q | phenomena. For bond percolation, the evidence of such a crossover of 

O ■ dimensionality would require an improvement of the GM formula to 

reach a relative error of typically 0.2%, while it is currently at 0.9% 
^ ; for hypercubes. 
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1 Introduction 



The discovery of renormalization group technics by Wilson in the early seven- 
ties has allowed the powerful elucidation of the mystery of critical phenomena 
PP. It is based on the existence of relevant variables, irrelevant variables and 
universality classes. Accordingly all parameters are classified as universal 
quantities and non universal quantities. For continuous phase transitions 
the critical exponents are universal while critical temperatures are not. 

In this framework dimension play a key role to categorize the effects of 
fluctuations. At very low dimensions, fluctuations are too strong and prevent 
any long range order to occur. The limit from which it does not happen is 
called the lower critical dimension d\. Only for d > d\ can long range order 
sustain fluctuations. On the other extreme, there exists some dimension d c 
called the upper critical dimension d c , beyond which fluctuations are averaged 
out and do not influence the critical properties. For d > d c there exists only 
one class of universality: the mean field one. Therefore it is in the range 
di < d < d c that fluctuations are instrumental to determine the critical 
properties. 

In parallel, percolation is a geometric phenomena with no temperature. 
However it was shown to be indeed identical to usual critical phenomena 
with d c = 6. Therefore its critical exponents are universal quantities while 
percolation thresholds p c are not. Accordingly the value of p c must be cal- 
culated for each system and varies from one geometry to another. However, 
at odd with this universality principle, a lot of efforts have been devoted to 
the finding of formulas for percolation threshold since about half a century. 
Several formulas have been proposed, which involve only two parameters : 
the dimension d and the coordination number q [21 El EJ El El El • The limits 
of such a choice has been discussed by J. C. Wierner et al [El E], who also 
considered that the most accurate of such formulas are the Galam-Mauger 
(GM) laws Indeed the high degree of accuracy of the GM law predictions 
hints for the existence of an underlying universality principle for percolation 
thresholds. 

In this work, we report on a possible crossover of a non universal quan- 
tity at the upper critical dimensionality in the field of percolation. Using 
the GM law and a series of recent numerical estimates for hypercube per- 
colation thresholds ( 6 < d < 13), site percolation thresholds are found to 
undergo a drastic change of behavior at the percolation upper critical di- 
mension d c = 6. This result is reminiscent of the crossover undergone by 
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universal quantities in critical phenomena. At contrast, nothing similar is 
evidenced at any dimension for the bond percolation. 

Considering only the case d > 2 of interest in the present work, the GM 
law is split in two. One applies to 3 < d < d c , and can be written: 

p c = P oi(d-l)(q-l)]~ a d b 3<d<d c . (1) 

The site percolation is approximated by Eq. (1), with b = 0, po = 
1.2868, a = 0.6160. As for the bond percolation threshold, po = 0.7541, b = 
a = 0.9346. This law corresponds to the so-called second class (the first one 
being related to d — 2 only), and will be called for this reason the GM2 law. 
Eq. (1) cannot extend up to d — > oo. Among several reasons outlined in [H], 
one comes from the fact that, in this limit, the percolation threshold should 
reduce to that of the Cayley tree for both sites and bonds. In other words, 
one must recover the Bethe asymptotic limit: 

Pc=pf = (9 -I)" 1 rf-^oo. (2) 

Eq. (2) is violated by the GM2 law. This drawback of the GM2 law was 
the main motivation for introducing another law associated to a third class, 
which applies at high dimension, and has the proper Bethe asymptotic limit. 
This is the asymptotic GM3 law: 

p c = 2 a - l [(d - l)(q - l)]- a d 2a -\ d » d c (3) 

with a = 0.08800 for sites, and a = 0.3685 for bonds. 

Despite the fact that the GM2 law is not exact, its accuracy is sufficient 
to materialize the dimensional dependence of the percolation threshold for 
the bcc, fee, hypercubic lattices up to d — 6. Due to the lack of data for the 
percolation thresholds in larger dimension, however, it was not possible so 
far to explore the existence of a critical crossover dimension d c above which 
the percolation threshold would follow a formula different from the GM2 law 
(approximated by he GM3 law). Recent Monte Carlo estimates for site and 
bond percolation thresholds with negligible standard deviations in simple 
hypercube lattices from d = 6 up to d = 13 ^U] now makes this investigation 
possible. 

It is the purpose of this work to make the comparison of the data including 
these Monte Carlo results with the predictions of Eq. (1) and (3). For site 
percolation threshold, the crossover between the two laws is clearly evidenced, 



3 



at dimension d c = 6. For bond percolation thresholds, however, no sizeable 
deviations from the GM2 law as defined by Eq. (1) is detected, up to the 
highest dimension d = 13 investigated. Such a crossover for bonds, if it exists, 
cannot be detected, since the GM2 and GM3 laws do not depart significantly 
from each other in the range 7 < d < 13. 

2 Analysis 

The numerical estimates pf , pf of the percolation thresholds for sites and 
bonds, respectively, are reported in Table 1, together with the results of the 
GM2 and GM3 laws. For sc lattices at d=5 and d=6, the data of refs. HUH2] 
have been substituted by those of ref . > since they are more accurate; we 
shall return to this point later on. 

2.1 site percolation 

As b = for site percolation, the GM2 law is best illustrated in a log-log plot 
of pf versus (d— l)(g — 1) in which case it is a straight line. This is illustrated 
in Fg. (1). We have reported on the same plot the numerical results of pf 
taken from El Il2] . For comparison, we have also reported (crosses) 
the values of pf predicted by the asymptotic GM3 law. Note according the 
GM3 law, (d — l)(q — 1) is not the pertinent variable, hence a random-like 
distribution of the crosses which cannot be connected to generate a curve. 
Indeed, since the pf 's depend on both d and q, only the crosses corresponding 
to lattices with the same topology, defined by the relation linking d and q 
can be connected. 

In practice, it means that the crosses corresponding to all the hypercubes 
(sc) from d = 3 up to d = 13 do belong to a single curve, since the same 
relation d = 2q holds true for all these lattices. The crosses corresponding 
to fee lattices should also belong to another curve, but the fee percolation 
threshold for site is known for d=3, 4 and 5 only, and three data points are 
not sufficient to materialize a curve. All the pf 's up to d=6 line up on the 
GM2 law (within the uncertainty limit above mentioned) as stated in ref. 

The new data for the sc lattices at higher dimensions, however, give 
evidence of a deviation of the pf 's from the linear GM law which increases 
with d, illustrated in Fig. (1). An equivalent formulation is to note a negative 
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lattice 


Pc(GM2) 


Pc(GM3) 


rrr 


Pc(GM2) 


Pc(GM3) 


TT7 \ 

P c \nu) 


Kagome 


0.65400 


0.59264 


0.65270 


0.51620 


0.38885 


0.52440 


Diamond 


[jAZb/O 


n /iQoo/i 
U.4ooz4 


U.4oUUU 


n on a c a 
U.oy4o4 


n o/ino/i 
U.z4yo4 




sc (d=3) 


0.31154 


0.27957 


0.31160 


0.24476 


0.20697 


0.24880 


bcc (d=3) 


0.25322 


0.20792 


0.24600 


0.17872 


0.18284 


0.18030 


fee (d=3) 


0.19168 


0.13969 


0.19800 


0.11714 


0.15479 


0.11900 


sc (d=4) 


0.19725 


0.18109 


0.19700 


0.16009 


0.14599 


0.16010 


fee (d=4) 


0.094794 


0.063570 


0.098000 








sc (d=5) 


0.14152 


0.13352 


0.14100 


0.11917 


0.11287 


0.11820 


fee (d=5) 


0.057352 


0.036740 


0.054000 








sc (d=6) 


0.10901 


0.10562 


0.10902 


0.095092 


0.092028 


0.094202 


sc (d=7) 


0.087898 


0.087319 


0.088951 


0.079233 


0.077697 


0.078675 


sc (d=8) 


0.073191 


0.074401 


0.075210 


0.067991 


0.067232 


0.067708 


sc (d=9) 


0.062410 


0.064802 


0.065210 


0.059602 


0.059255 


0.059496 


sc (d=10) 


0.054198 


0.057390 


0.057593 


0.053097 


0.052971 


0.053093 


sc (d=ll) 


0.047756 


0.051497 


0.051590 


0.047904 


0.047893 


0.047795 


sc (d=12) 


0.042578 


0.046699 


0.046731 


0.043659 


0.043704 


0.043724 


sc (d=13) 


0.038336 


0.042718 


0.042715 


0.040124 


0.040190 


0.040188 



Table 1: Numerical estimates for the percolation thresholds pf (nu) , together 
with the results of the GM2 and GM3 laws. For sc lattices at d=5 and d=6, 
the data of refs. [HJH2] have been substituted by those of ref. jTUj, since they 
are more accurate. 
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curvature of the 'numerical' Pc(d) curve for sc hypercube lattices at d > 6 
to match the GM3 law. To quantify this effect, we have plotted in Fig. (2) 
the relative difference Apf/pf = [p^{nu) — pf (GM2)]/pf {nu) as a function 
of d for the hypercubes (since data are available only for these lattices at 
high dimensions). p'^(nu) is the numerical percolation threshold JU] and 
pf (GM2) the prediction of the GM2 law. 

We have reported in [U] that \Ap^\ can reach 0.008 for some lattices in 
d < 6. This measures the accuracy of the GM2 law when applied to any 
Bravais lattice sc, be, fee in 3 < d < 6. However, regarding the simple 
cubic and hypercubes only, the accuracy is much better. In particular, in 
our prior work (see Table 1 in ref. 6 ), |Apf | for sc lattices was within 
5xl0~ 4 at all dimensions d < 7, except at d = 6 where an outstanding 
deviation |Apf | = 0.002 was pointed out between the numerical estimate 
0.107 available at that time, and 0.109 predicted by the GM2 law. The 
new numerical calculations ^U] have corrected the estimate of pf at d — 6, 
raising from 0.107 to 0.109017, now in agreement with the prediction of 
the GM2 law. With the new estimates of ref. ^U] which are ~ 30 times more 
precise than the previous ones, we find |Apf|/pf < 0.4% in the whole range 
3 < d < 6 That is why we consider as significant a departure from the 
GM2 law with |Apf |/pf > 0.4% in Fig. (2). 

This plot then provides evidence for a crossover at d c = 6: at d < 6, the 
GM2 law applies; at d > 6, this is no longer the case. The systematic quasi 
linear increase of Ap^ jp s c as a function of d which extrapolates to zero at 
d = 6 corroborates this value for the crossover dimension. Note that d c is also 
known to be the upper marginal dimension where critical exponents for the 
percolation transition reach mean-field values. However, in phase transition 
theory, only universal quantities are supposed to undergo a crossover at d c . 
This site percolation transition gives an outstanding example where a non 
universal quantity like pf also undergoes a crossover at d c . 

Instead of choosing the GM2 law as the reference, we can also choose 
the GM3 law, and investigate how the the percolation thresholds approach 
this asymptotic law. The GM3 law is best illustrated in the log-log plot 
of 2dp s c versus x = 2d 2 /[(d — l)(q — 1)] in Fig. (3) since it reduces to a 
straight line. For comparison, the data have been also reported, and the 
crosses now correspond to the values of pf as predicted by the GM2 law. 
As x is not the pertinent variable according to the GM2 law, once again, 
only the crosses corresponding to all the hypercubes (sc) from d = 3 up to 
d — 13 do belong to a single curve. As we can see on fig. 2, this curve has 
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a negative curvature and crosses the straight line corresponding to the GM3 
law. pf(nu) as a function of d then shifts from the GM2 law at d c = 6, to 
approach the GM3(law) assumed to be its asymptote in the GM model. 

To illustrate this behavior, we have reported in Fig. (4) Apf/pf = 
\pf(nu) — pf (GM3)]/pf as a function of d for hypercubes in high dimen- 
sions. From this Figure, it can be seen that this asymptotic limit is indeed 
reached at d — 13. It is then important to note that the crossover at d c = 6 
does not mean an abrupt shift from the GM2 law to the GM3 law. Instead, 
it is a crossover to another law which is missing here. More important, it 
indicates this law pf (nu) as a function of d is not embedded in the GM for- 
mula, and accepts the GM3 law only as an asymptote in the large d limit, 
eventually reached (within the error bars) at d > 13. Actually, the system- 
atic and increasing deviation of pf from the GM3 law as d decreases from 
d — 13 can be viewed as a pretransitional effect upon approaching the upper 
critical dimension from below, beyond the scope of the GM3 law. 

2.2 Bond percolation 

Let us now investigate the situation for bonds. The GM2 law in this case 
is illustrated in a log- log plot of pf as a function of (d — l)(q — \)/d (Fig. 
(5)). No deviation from the GM2 law can be evidenced for any lattice up 
to the highest dimension d — 13 investigated. To be more specific, we have 
reported in Fig. (6) the differences Apf /pf with Apf = pf (GM2) —pf(nu) 
and Apf? = pf (rai) — pf (GM3). Note the sign inversion in the definition of 
Apf to have this quantity positive, which means that the numerical estimates 
are in between the estimates of the GM2 and the GM3 laws. Clearly, the 
exact (numerical) results match the GM2 law at d < 6, as expected. The 
relative deviation of the GM2 law with respect to the numerical result at 
d < 6 is in the range Apf/pf < 0.9%, so that we can take as significant a 
departure from the GM2 law Apf/pf > 0.9%. 

However, the relative deviation between numerical estimates and the GM2 
law is smaller in all dimensions 7 < d < 13. Actually, the relative deviation of 
the numerical estimates does not exceed 0.4% in the whole range 9 < d < 13, 
not only with respect to the GM2 law, but also with respect to the GM3 
law, so that no significant difference between the two laws can be detected. 
Therefore, at contrast with the situation met in the site percolation problem, 
the GM3 law is not accurate enough to detect a crossover at d — 6. 

To explain why it is not possible to distinguish between the GM2 and 
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the GM3 laws in the range 7 < d < 13, we note that the leading term in 
a q — 1 expansion of the GM2 law at large q is in (q — l) b ~ 2a . For bonds, 
b — 2a = —0.94, very close to the exponent -1 of the leading term in the GM3 
law. For sites, however, b — 2a = —1.23 which is markedly different from -1, 
so that a clear distinction between the GM2 law and the GM3 law can be 
made even at the scale of a short dimension interval for sites, while this is 
impossible for bonds. 

To be more specific, we note that the deviation of both the GM2 and 
GM3 laws with respect to the numerical results, i.e. \p^(GM2) — pf (nu)\ 
and (GM3) — pf (nu)\, is only on the fourth digit in dimension d > 7, 
which amounts to a relative error within 0.4%, while, as, we have stated 
above, a deviation with respect to the GM laws can be regarded as significant 
only if it exceeds 0.9% for bonds. Therefore, the evidence of a crossover of 
dimensionality at d c — 6 for bonds requires an improvement of the GM law 
and its substitution by a new formula which would improve the accuracy 
typically a factor 3 

3 Discussion 

Van der Marck has shown that, if there is to be an exact universal 
formula for percolation thresholds, it must be based on more information 
than d and q only. In particular, the body centered cubic lattice and the 
stacked triangular lattice both realize d = 3, q = 8, but they have different 
p c 's. This simple consideration is sufficient to show that Eqs. (1, 3) cannot 
be exact. Nevertheless, the deviation with respect to the exact or the almost 
exact (deduced from the best numerical estimates) thresholds of sc, bcc and 
fee lattices is so small that the GM laws have the operational power to address 
the problem of the crossover dimensionality. 

On this basis, the present analysis of the percolation thresholds for sc (hy- 
percubes) hints at the possible existence of a crossover at d c — 6 for the site 
percolation threshold, although it is a non-universal parameter. At higher 
dimensions, departs from the GM2 law to approach the GM3 asymptotic 
law which is reached at d — 13. For bonds no similar behavior at d c — 6 is 
detected. 

Because of the relation d = 2q for hypercubes, the study of hypercubes 
alone does not allow us to distinguish between the variables d and q to know 
which one is pertinent for the possible crossover. However, some indication 
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can be extracted from the data in Fig. (1). Let us consider in particular the 
fee lattice, which has a coordinance q much larger than the hypercube at the 
same dimension. For instance, q = 24 and 40 for fee in dimensions d = 4 and 
d — 5, respectively. Those are the coordination numbers of hypercubes in 
dimensions 12 and 20, respectively. The fee percolation thresholds at d — 4, 5 
have been reported in Fig. (1), along with the values predicted by the GM2 
and GM3 laws. This Figure shows unambiguously that these two fee lattices 
do satisfy the GM2 law and not the GM3 law, despite the fact that their 
coordinance number are those of the hypercubes in dimensions d — 12 and 
d= 14. 

It corroborates that the pertinent variable responsible for the possible 
crossover behavior in the percolation thresholds is indeed the dimension d, 
while the coordination q does not play any significant role. It also justifies the 
term crossover of dimensionality used in this work. For universal variables, 
such as the critical exponents, the renormalization group technique provides 
us with an elaborate theory to understand the crossover of dimensionality, at 
the associated identification of universality classes. However such a behavior 
is not expected for non-universal quantities like percolation thresholds. The 
explanation for this unexpected result is thus a new challenge in the field of 
phase transition theory. 

We note that d c is the upper critical limit for percolation transition phe- 
nomena in general, not for the hypercubes only. The question then arises 
whether the possible crossover of dimensionality we have observed for hyper- 
cubes also applies to other lattices with a different topology. So far, accurate 
numerical estimates of percolation thresholds beyond d — 6 could be achieved 
on the hypercubes only. The lack of data prevents us for the moment from 
addressing this question. Nevertheless, we believe that the progress in the 
methods to compute percolation thresholds will make possible the simula- 
tion of percolation of other systems more complex with a larger coordinate 
number q in the near future. 
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Figure 1: Site percolation threshold as a function of the variable (d-l)(q-l) 
pertinent to the GM2 law, in decimal logarithm. The numerical estimates (•) 
are from refs. |lUtlTT| ll2j. The crosses are predictions of the GM3 asymptotic 
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Figure 2: Relative difference [pf (rax) — pf (GM2)]/pf (n«) between the numer- 
ical data pf (rax) and the the GM2 law pf (GM2) (full squares) clS cl function 
of the dimension ci of the hypercubes. Lines are guides for the eyes. 
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Figure 3: Site percolation threshold as a function of the variable 
2d 2 /[(d— l)(q — 1)] pertinent to the GM3 law, in decimal logarithm. The 
numerical estimates (•) are from refs. |lUtlTH ll2j. The crosses are predictions 
of the GM2 asymptotic law, the solid line materializes the GM3 law. 
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Figure 4: Relative difference [pf(mx) — pf (GM3)]/pf {nu) between the nu- 
merical data pf(m/)and the the GM3 law p^(GM3) as a function of the 
dimension d of the hypercubes. 
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Figure 5: Bond percolation threshold as a function of the variable 
[d — l)(g — pertinent to the GM2 law, in decimal logarithm. The nu- 
merical estimates (•) are from refs. |lut ITT| IT2"]. The solid line materializes 
the GM2 law. 
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Figure 6: Relative difference Apf/pf (mi) with Apf = pf(GM2) —p^(nu) 
(crosses) and Apf = pf(nu) — pf(GM3) ( full squares) as a function of the 
dimension d of the hypercubes. 
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